Observation of the Naive-T-odd Sivers Effect in Deep-Inelastic Scattering 
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Azimuthal single-spin asymmetries of lepto-produced pions and charged kaons were measured 
on a transversely polarized hydrogen target. Evidence for a naive- T-odd, transverse-momentum- 
dependent parton distribution function is deduced from non- vanishing Sivers effects for ty + , 7r°, and 
K , as well as in the difference of the 7r + and n~ cross sections. 
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The ongoing experimental effort in spin-dependent 
high-energy scattering and attendant theoretical work 
continue to indicate that the spins of the quarks and glu- 
ons are not sufficient to explain the nucleon spin [J] . The 
investigation of the only remaining contribution, that of 
orbital angular momentum of the constituents, is clearly 
essential. Transverse-momentum-dependent parton dis- 
tribution functions are recognized as a tool to study spin- 
orbit correlations, hence providing experimental observ- 
ables for studying orbital angular momentum. One par- 
ticular example is the Sivers function f^ T describing 
the correlation between the momentum direction of the 
struck quark and the spin of its parent nucleon. This 
correlation is commonly defined as the Sivers effect. A 
non- vanishing f^p contributes to, e.g., single-spin asym- 
metries (SSAs) in semi-inclusive deep-inelastic scattering 
(DIS) off transversely polarized protons, ep^ — > e'hX, 
where h is a hadron detected in coincidence with the 
scattered lepton e'. 

For a long time, transverse SSAs had been assumed 
to be negligible in hard scattering processes. They are 
odd under naive time reversal, i.e., time reversal of three- 
momenta and angular momenta, and thus require inter- 
ference of amplitudes with different hclicities and phases. 
In QED and perturbativc QCD, these ingredients are 
suppressed 0, ij. Therefore, in semi-inclusive DIS they 
must be ascribed to the non-perturbative parts in the 
cross section, i.e., to specific parton distribution and 
fragmentation functions, commonly categorized as being 
naive- T-odd. The idea of a naive- T-odd quark distribu- 
tion function goes back to an interpretation Q of large 
left-right asymmetries observed in pion production in the 
collision of unpolarized with transversely polarized nuclc- 
ons 0. It was argued that such asymmetries could be 
attributed to a left-right asymmetry in the distribution 
of unpolarized quarks in transversely polarized nucleons, 
i.e., an asymmetry that exists before the pion is formed 
in the fragmentation process, and that does not vanish at 
high energies. A decade after an initial proof [f| that this 
distribution function, now termed the Sivers function, 
must vanish because of time-reversal invariance of QCD, 
it was realized through the pioneering work in Ref. Q 
and subsequently in Refs. [a S El that this proof ap- 
plies only to transversc-momcntum-integrated distribu- 
tion functions. A gauge link, previously neglected in 
the definition of gauge-invariant distribution functions, 
invalidates the original proof for the case of transverse- 
momcntum-dependent distribution functions. The gauge 
link provides the phase for the interference (required for 
naive-T-oddncss) , and can be interpreted as an interac- 
tion of the struck quark with the color field of the target 
remnant [ill ]. 

The inclusion of the gauge link has profound conse- 
quences on factorization proofs and on the concept of 
universality, which are of fundamental relevance for high- 
energy hadronic physics. A direct QCD prediction is a 



Sivers effect in the Drell-Yan process that has the oppo- 
site sign compared to the one in semi- inclusive DIS [8[. 
For hadron production in pro ton-proton collisions the sit- 
uation is more intricate [12j, leading to a violation of 
standard factorization and universality, even for the case 
of unpolarized collisions [13j . Therefore, the study of the 
Sivers effect in semi-inclusive DIS and other processes 
is of utmost importance for our understanding of high- 
energy scattering involving hadrons. 

The Sivers effect has been related to the orbital motion 
of quarks inside a transversely polarized nucleon since the 
seminal work in Ref. ■ In the calculation of Ref. [3] , it 
became clear that orbital angular momentum of quarks 
is needed for a non-vanishing Sivers effect as it arises 
through overlap integrals of wave-function components 
with different orbital angular momenta. However, no 
quantitative relation has yet been found between f^p and 
the orbital angular momentum of quarks. One faces a 
similar quandary with the anomalous magnetic moment 
k of the nucleon: it also requires wave function compo- 
nents with non-vanishing quark orbital angular momen- 
tum without constraining the net orbital angular momen- 
tum 14 1 . Indeed, f^ T involves overlap integrals between 
the same wave function components that also appear in 
the expressions for k as well as for the total angular mo- 
mentum in the Ji relation [l5| for the nucleon-spin de- 
composition 0, EH. 

An interesting link between k and f^p was suggested 
in Ref. [l6j]: the sign of the quark-flavor contribution to 
k determines the sign of f^p for that quark flavor. If the 
final-state interactions are attractive, as one would as- 
sume for the confining color force, a positive flavor con- 
tribution to k leads to a negative /fj.. (The sign and 
angle definitions follow the Trento Conventions |l7|.) 

In semi-inclusive DIS, f^p leads to SSAs in the dis- 
tribution of hadrons in the azimuthal angle about the 
virtual-photon direction. In general, azimuthal SSAs 
provide important information not only about the Sivers 
function but also about other distribution and fragmen- 
tation functions. For example, transversity [18[, de- 
scribing the distribution of transversely polarized quarks 
in transversely polarized nucleons, combined with the 
naive- T-odd Collins fragmentation function Q , also leads 
to SSAs. The keys to extracting different combinations 
of the various distribution and fragmentation functions 
are their different dependences on the two azimuthal an- 
gles <j> and <ps of the hadron momentum Ph and of the 
transverse component St of the target-proton spin, re- 
spectively, about the virtual-photon direction (cf. [171]). 
The Sivers effect manifests itself as a sin(0 — <j>s ) modu- 
lation in the azimuthal distribution (ljj . 

In this Letter clear evidence for a non- vanishing Sivers 
function is reported. The sin((f> — <pg) modulations in 
semi-inclusive DIS are measured for pions and charged 
kaons, as well as in the difference between the ir + and 
7r~ cross sections, providing sensitivity to f^p for both 
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valence and sea quarks. 

The data reported here were recorded during the 2002- 
2005 running period of the Hermes experiment using 
a transversely nuclear-polarized hydrogen gas target in- 
ternal to the 27.6 GeV Hera lepton (e + or e~) storage 
ring at Desy. The open-ended target cell was fed by an 
atomic-beam source (2(| based on Stern-Gerlach sepa- 
ration combined with radio-frequency transitions of hy- 
pcrfine states. The nuclear spin direction was flipped 
at 1-3 min time intervals, while both nuclear polariza- 
tion and the atomic fraction inside the target cell were 
continuously measured [2l|. The average magnitude of 
the proton-polarization component perpendicular to the 
lepton-beam direction was 0.725 ± 0.053. 

Scattered leptons and coincident hadrons were de- 
tected by the Hermes spectrometer 22j. Leptons were 
identified with an efficiency exceeding 98% and a hadron 
contamination of less than 1%. Charged hadrons with 
momentum 2 GeV < \Ph\ < 15 GeV were identified 
using a dual-radiator ring-imaging Cerenkov detector 
(RICH) [23j ]. For this a hadron-identification algorithm 
was employed that takes into account the topology of 
the whole event, in contrast to the track- level algorithm 
in previous analyses Events were selected sub- 

ject to the requirements Q 2 > lGeV 2 , W 2 > 10 GeV 2 , 
O.K y < 0.95, and 0.023 < x < 0.4, where Q 2 = -q 2 = 
-(k - k'f, W 2 = (P + q) 2 , y = (P ■ q)/(P ■ k), and 
x = Q 2 /{2P ■ q). Here, P, k, and k' represent the four- 
momenta of the target proton, the incident lepton, and 
the scattered lepton, respectively. Coincident hadrons 
were accepted if 0.2 < z < 0.7, where z = (P ■ Ph)/ (P -q). 

The cross section for semi-inclusive production of 
hadrons using an unpolarized lepton beam on a trans- 
versely polarized target can be written as (lil . 25 . 2r| 



a{4>, 4>s) = cr uu {l + 2(cos 4>) uu cos + 2(cos 2(j)) uu cos 20 
+ |S T | [2(sin(0-^)) UT sin(0 -<t> s ) + .. ■}}, (1) 

where five sine modulations contribute to the 
polarization-dependent part, but, for convenience, 
only the sin(0 — (j)s) modulation (the Sivers term), is 
written out explicitly. Here, the subscript UT denotes 
unpolarized beam and transverse target polarization 
(with respect to the virtual-photon direction), while a ul7 
represents the ^-independent part of the polarization- 
independent cross section. The sin(0 — <j>s) amplitude 
can be interpreted in the quark-parton model as 



1!: 



2(sm((j)-(j)s)) ul 



Z q e 2 q f!(x,p 2 T ) ®D\{z,K 2 ) 



(2) 

where the sums run over the quark flavors, the e q are the 
quark charges, and j\ and D\ are the spin-independent 
quark distribution and fragmentation functions, respec- 
tively. The symbol ® (®w) represents a (weighted) con- 
volution integral over intrinsic and fragmentation trans- 
verse momenta p T and Kt, respectively 



The amplitudes of the five sine modulations in Eq. (p} 
were extracted simultaneously to avoid cross contamina- 
tion. For this a maximum- likelihood fit was used [27j |. 
with the data alternately binned in x, z, and Phi. = 
\Ph — \ q $ 1) but unbinned in and <ps- A sixth term, 
arising from the small but non- vanishing target-spin com- 
ponent that is longitudinal to the virtual-photon direc- 
tion when the targ et is polarized perpendicular to the 
beam direction [2a |. was also included in the fit. 

A scale uncertainty of 7.3% on the extracted Sivers 
amplitudes arises from the accuracy of the target- 
polarization determination. Inclusion in the fit of es- 
timates [2^] for the cos and cos 20 amplitudes of the 
unpolarized cross section had negligible effects on the 
amplitudes extracted. Possible contributions [28| to the 
amplitudes from the non-vanishing longitudinal target- 
spin component were estimated based on measurements 
of SSAs on longitudinally polarized protons 3(J 31] and 
included in the systematic uncertainty. Effects from the 
hadron identification using the RICH, the geometric ac- 
ceptance, smearing due to detector resolution, and radia- 
tive effects are not corrected for in the data. Rather, the 
size of all these effects was estimated using a simulation 
tuned to the data, which involved a fully differential poly- 
nomial fit to the measured azimuthal amplitudes [321 ] . 
The result was included in the systematic uncertainty 
and constitutes the largest contribution. 

Based on a Pythia6 Monte Carlo simulation [33| 
tuned to Hermes data, the fraction of charged pions 
(kaons) stemming from the decay of exclusive vector- 
meson channels was estimated to be about 6-7% (2-3%). 
Among the contributions of all the vector mesons to the 
pion samples, that of the p° is dominant. A different 
observable, for which the contributions from exclusive p° 
mesons cancels, is the pion- difference asymmetry 



Al, T * (0,0s) 



1 Kt-<)-(^ 
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S T \ (^;-a^) + (a^-^T) 



(3) 



the SSA in the difference in the tt + and ir~ cross sec- 
tions for opposite target-spin states In addition, 
this asymmetry helps to isolate the valence-quark Sivers 
functions: under some assumptions, such as charge- 
conjugation and isospin symmetry among pion fragmen- 
tation functions, one can deduce from Eq. ^ that this 
SSA stems mainly from the difference (f^ d " — 4/^"") 
in the Sivers functions for valence down and up quarks. 

The resulting Sivers amplitudes for pions, charged 
kaons, and for the pion-diffcrence asymmetry are shown 
in Fig. [1] as functions of x, z, or Phx- They are positive 
and increase with increasing z, except for 7r~, for which 
they are consistent with zero. In the case of tt + , K + , 
and the pion-diffcrence asymmetry, the data suggest a 
saturation of the amplitudes for Phi. > 0.4 GeV and are 
consistent with the predicted linear decrease in the limit 
of Ph± going to zero. 
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FIG. 1: Sivers amplitudes for pions, charged kaons, and the 
pion-difference asymmetry (as denoted in the panels) as func- 
tions of x, z, or Ph±- The systematic uncertainty is given as a 
band at the bottom of each panel. In addition there is a 7.3% 
scale uncertainty from the target-polarization measurement. 



In order to further examine the influence of exclusive 
vector-meson decay and other possible ^--suppressed 
contributions, several studies were performed. Raising 
the lower limit of Q 2 to 4 GeV 2 eliminates a large part 
of the vector-meson contribution. Because of strong cor- 
relations between x and Q 2 in the data, this is presented 
only for the z and Ph± dependences. No influence of the 
vector-meson fraction on the asymmetries is visible as 
shown in Fig. [21 For the x dependence shown in Fig. [3J 
each bin was divided into two Q 2 regions below and above 
the corresponding average Q 2 ({Q 2 (xi))) for that x bin. 
While the averages of the kinematics integrated over in 
those x bins do not differ significantly, the (Q 2 ) values 
for the two Q 2 ranges change by a factor of about 1.7. 
The asymmetries do not change by as much as would 
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FIG. 2: Sivers amplitudes for ir (left) and K (right) as 
functions of z or Ph± , compared for two different ranges in Q 2 
(high-Q 2 points are slightly shifted horizontally). The corre- 
sponding fraction of pions and kaons stemming from exclusive 
vector mesons, extracted from a Monte Carlo simulation, is 
provided in the bottom panels. 



c 0.1 

'jn 

CM 





<iT 10 
> 

91 
(3 



■ Q 2 < <Q 2 ( Xi )> 

o Q 2 > <Q 2 (x.)> 



K + 



10 



10 



FIG. 3: Sivers amplitudes for ty + (left) and K + (right) as 
functions of x. The Q 2 range for each bin was divided into 
the two regions above and below (Q 2 (xi)) of that bin. In the 
bottom the average Q 2 values are given for the two Q 2 ranges. 



have been expected for a sizable ^--suppressed contri- 
bution, e.g., the one from longitudinal photons to the 
spin- (in) dependent cross section. However, while the tt + 
asymmetries for the two Q 2 regions are fully consistent, 
there is a hint of systematically smaller K + asymmetries 
in the large-Q 2 region. 

An interesting facet of the data is the difference in the 
7r + and K + amplitudes shown in Fig. [4] On the basis 
of u-quark dominance, i.e., the dominant contribution to 
7r + and K + production from scattering off u-quarks, one 
might naively expect that the 7r + and K + amplitudes 
should be similar. The difference in the ir + and K + 
amplitudes may thus point to a significant role of other 
quark flavors, e.g., sea quarks. Strictly speaking, even in 
the case of scattering solely off u-quarks, the fragmenta- 
tion function D±, contained in both the numerator and 
denominator in Eq. (J2]), does not cancel in general as it 
appears in convolution integrals. This can lead not only 
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FIG. 4: Difference of Sivers amplitudes for K + and 7r + as 
functions of x for all Q 2 (left), and separated into "low-" and 
"high-Q 2 " regions as done for Fig. 



to additional z-dependences, but also to a difference in 
size of the Sivers amplitude for 7t + and K + . Higher- twist 
effects in kaon production might also contribute to the 
difference observed: in the low-Q 2 region, where higher- 
twist should be more pronounced, the ir + and K + am- 
plitudes disagree at the confidence level of at least 90%, 
based on a Student's t-test, while being statistically con- 
sistent in the high-Q 2 region. 

As scattering off u-quarks dominates in these data, the 
positive Sivers amplitudes for tt + and K ± suggest a large 
and negative Sivers function for u-quarks. This is sup- 
ported by the positive amplitudes of the difference asym- 
metry, which is dominated by the contribution from va- 
lence u-quarks. The vanishing amplitudes for 7r~ require 
cancelation effects, e.g., from a d-quark Sivers function 
opposite in sign to the it-quark Sivers function. In combi- 
nation with deuteron data from the COMPASS collabora- 
tion [34j , a large positive <i-quark Sivers function can be 
deduced [H( . These fits have yet to be updated with the 
final results presented here, as well as with preliminary 
proton data from Compass (36|. 

In summary, non-zero Sivers amplitudes in semi- 
inclusive DIS were measured for production of 7r + , 7r°, 
and K^, as well as for the pion-difference asymme- 
try. They can be explained by the non-vanishing naive- 
T-odd, transverse-momentum-dependent Sivers distribu- 
tion function. This function also plays an important 
role in transverse single-spin asymmetries in pp collisions, 
and is linked to orbital angular momentum of quarks in- 
side the nuclcon. Although no quantitative conclusion 
about their orbital angular momentum can be inferred, 
the Sivers function provides important constraints on the 
nucleon wave function and thus indirectly on the total 
quark orbital an gula r momentum. For instance, in the 
approach of Ref. , the measured positive Sivers asym- 
metries for 7r + and K + mesons correspond to a positive 
contribution of it-quarks to the orbital angular momen- 
tum, under the assumption that the production of 7T" 1 " 
and K + mesons is dominated by scattering off it-quarks. 
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